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Big Bang Nucleosynthesis and Active-Sterile Neutrino Mixing: 

Evidence for Maximal v T Mixing in Super Kamiokande? 

Xiangdong Shi and George M. Fuller 

Department of Physics, University of California, San Diego, La Jolla, California 92093-0319 
We discuss Big Bang Nucleosynthesis constraints on maximal «-> v s 

mixing. Vacuum f„ <-> v s oscillation has been proposed as one possible ex- 
planation of the Super Kamiokande atmospheric neutrino data. Based on the 
most recent primordial abundance measurements, we find that the effective 



O 

Q . number of neutrino species for Big Bang Nucleosynthesis (BBN) is N v < 3.3. 



Assuming that all three active neutrinos are light (with masses <C 1 MeV), 



^NJ ■ we examine BBN constraints on v„ <-> v s mixing in two scenarios: (1) a neg- 

. ligible lepton asymmetry (the standard picture); (2) the presence of a large 

. lepton asymmetry which has resulted from an amplification by v T «-> v s > mix- 

ing {y s i being v s or another sterile neutrino species). The latter scenario 

ON 

. has been proposed recently to reconcile the BBN constraints and large-angle 

q . <-> v s mixing. We find that the large-angle <-> v s mixing in the first see- 

nario, which would yield N v ~ 4, is ruled out as an explanation of the Super 
Kamiokande data. It is conceivably possible for the i/„ <-> u s solution to evade 
BBN bounds in the second scenario, but only if 200 eV 2 < m 2 r — m 2 , ; < 10 
eV 2 is satisfied, and if v T decays non-radiatively with a lifetime < 10 3 years. 
This mass-squared difference implies 15 eV< m Vr < 100 eV if v s i is much 
lighter than v T . We conclude that maximal (or near maximal) <->• v T mix- 
ing is a more likely explanation of the Super Kamiokande data. 

PACS numbers: 14.60.Pq; 14.60.St; 26.35.+c; 95.30.Cq 
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I. INTRODUCTION 



The recent Super Kamiokande (SuperK) data on atmospheric neutrinos show strong 
evidence for neutrino oscillation in a channel which involves muon neutrinos and another 
neutrino species, and with vacuum mass-squared difference 3 x 10~ 4 eV 2 < 5m 2 < 7 x 10 -3 
eV 2 and vacuum mixing angle satisfying 0.8 < sin 2 29 < 1 [|T],|2|]. Measurements of the 
muon-like events and limits on the electron-like events preclude mixing between and v e 
as a significant channel for reduction M. A <-> u e mixing with these parameters 
would also yield a uniform suppression of the solar neutrino flux by roughly a factor ~ 2 
across the entire solar neutrino energy spectrum. However, such an energy independent 
suppression contradicts what is implied from the combined solar neutrino data [[|. This 
leaves two maximal or near maximal vacuum neutrino oscillation channels as possible avenues 
for solution: <-> v T mixing; or v s mixing, where v s is a sterile neutrino (e.g., a singlet 

under the SU(3)x SU(2)x U(l) gauge symmetry of the Standard Model). 

However, it has been pointed out by many groups that z/ M «-> v s mixing with the 
aforementioned parameters would violate the constraints from Big Bang Nucleosynthesis 
(BBN) if all three active neutrinos are light compared to 1 MeV. This violation results 
from <-> v s mixing essentially bringing an extra v s into chemical equilibrium with the 
three known species at the BBN epoch. In turn, these extra degrees of freedom imply a 
larger expansion rate and so would yield a primordial 4 He abundance that is too high to 
accommodate observationally-determined abundance bounds 0-0. 

There has been some confusion on this result lately due to a larger systematic uncertainty 
in the measured primordial 4 He abundance as argued by some groups, and due to recent 
measurements of the primordial deuterium abundance. Furthermore, the BBN constraint 
has been obtained by assuming a negligible lepton number asymmetry (< 10~ 7 ), which has 
been shown to be an invalid assumption under circumstances where the lepton asymmetries 
in the neutrino sector can be amplified by active-sterile neutrino transformations |§|9]]. It 
has been argued that a v T mixing with a lighter sterile v s i may amplify an initially negligible 
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lepton asymmetry to a large level that either is sufficient to suppress a subsequent <-> v s 
mixing j5|, or can be partly converted to a positive chemical potential in the v e v e sector to 
lower the effective N v ||10|| . Both may conceivably lead to an evasion of the BBN bound on 
the <-» v s mixing parameters. 

It is traditional in BBN studies to parametrize both the 4 He yield and the expansion rate 
of the universe at the BBN epoch in terms of an effective number of relativistic neutrino 
flavors, N u . In the context of neutrino oscillation discussions this convention is potentially 
confusing and misleading, since efficient neutrino matter-enhanced transformation could 
result in both a larger neutrino energy density and a lower 4 He yield (by way of a v e v e 
asymmetry) than in the standard picture! Here we will also employ N v in the traditional 
manner to facilitate comparison with the results of previous work. However, we will attempt 
to point out what the true underlying picture is in every case. 

In this paper we intend to clarify the bound from BBN in light of the latest development 
in the primordial 4 He and D abundance measurements, and critique the possibility that the 
Vp, «-> v s mixing parameters required to explain the SuperK data might survive the BBN 
bound because of an amplified lepton asymmetry. We implicitly assume that all three active 
neutrinos are light, being relativistic at the BBN epoch. 

II. FORMALISM 

Here we briefly outline the formalism for active-sterile neutrino transformation used in 
this paper. A detailed treatment of active-sterile neutrino mixings in the BBN epoch can 
be found in ref. or 0. 

We adopt units in which h = c = k = 1. We denote by n« the number density of a 
particle % relative to its equilibrium value at a temperature T. These equilibrium values 
are 2((3)T 3 /n 2 for photons, 3((3)T 3 /2ir 2 for electrons/positrons, and 3C(3)T 3 /47r 2 for each 
neutrino species. We describe the u a <-> v s (a = e, fi or r) transformation channel at the 
BBN epoch with the following two differential equations: 
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(T(v a v a — ► ii))(nini - n Ua n p J; (1) 



In the equations, Pq denotes the total number density of the mixture of active and sterile 
neutrino species, Pq = n Ua +n Us , and P denotes the composition of the mixture. In particular, 
Pz — n v a — n u s - The other component of P, P^ = -P^x + PyYi is an indication of the phase 
coherence of the oscillation channel. Quantities (T) are thermally averaged reaction rates. 
The vector V represents the frequency and the axis of the oscillation in P-space, and the D- 
term represents the damping of P^ due to neutrino interactions. This latter term constantly 
acts to reduce a mixed neutrino state into flavor eigenstates of v a and u s . 
At the epoch of BBN, 

V x = —wn2ff, V y = 0, V z = -— C os29 + V^ + V^ (3) 

where 5m 2 = m 2 Ug — m 2 a and 9 are the usual vacuum mixing parameters of mass-squared 
difference and effective vacuum 2x2 mixing angle (5m 2 > if v s is heavier than u a ), and E 
is the energy of the neutrinos. is the contribution of the MSW matter-enhanced effect 



driven by asymmetries in the background plasma |TT|, and is approximately 



Vt « 0.13G F T 3 [8L /3 + 2(n„„ - n Da ) + ]T (n V0 - n^)] (4) 

where Gf is the Fermi constant, T is the plasma temperature, and Lq represents the con- 
tributions from the baryonic asymmetry and electron-positron asymmetry, ~ 10~ 9 . is 



the contribution due to finite temperature neutrino mass renormalization effects | TT| , and 
for a = fi or r, is 

V^-ll(n Ua +n Pa )G 2 F ET\ (5) 
Meanwhile, for a = fi or r, we also have 

D ~ 0.5G|T 4 £. (6) 



It has been shown that adopting the average neutrino energy E m 3.151T in Eq. and 
thermally averaging interaction rates give a fairly good description of active-sterile neutrino 



transformations in BBN []E|. 

The same equations ([!]) through ([5]) also apply to the u a v s mixing if we switch n Va 
(n V/3 ) and n 9a (n 9/3 ) in the expression of V£. Here we add an overbar on top of all variables 
associated with the u a <-> v s system to denote the corresponding variables in the v a <-> v s 
system. 

The initial condition of the universe is taken to be a pure v a v a system (i.e., Pq = Pq = 
P z = P z = 1 and P x = P x = P y = P y = 0) at an initial temperature 

5m 2 cos 29 



T init » 15 



1/6 

(7) 



eV 2 

This temperature is taken to be high enough that Vj dominates over Sm 2 /2E, guaranteeing 
that any neutrino transformation is severely suppressed. (Note however, that Ij n it has to be 
below the temperature of the Quark-Hardron phase transition, ~ 150 MeV ||13|| .) 

Barring the existence of a large initial lepton asymmetry generated by unknown processes 
at T > 150 MeV, we will assume a negligible initial lepton asymmetry L (e.g., at the 
same level as L ). Here a significant L can be generated only by active-sterile neutrino 
transformation during the BBN epoch. 

III. BBN BOUND WHEN LEPTON ASYMMETRY IS NEGLIGIBLE 

The primordial 4 He yield Y in standard BBN is principally a function of the expansion 
rate of the universe, as this determines the temperature of weak freeze-out and, hence, the 
overall neutron-to-proton ratio. Secondarily, the 4 He abundance also depends on the baryon 
density iVe (usually expressed as rj = N-q/N^), albeit only weakly. Therefore, given an 
7] inferred from the primordial deuterium abundance D/H (which depends sensitively on 
rj), a bound on N v in the standard BBN picture can be obtained based on the measured 
primordial 4 He abundance ||14||. The same bound can be used to constrain active-sterile 
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neutrino mixings, because such mixings can produce extra neutrino energy density. However, 
one should use caution here because these scenarios could also result in a larger number of 
u e 's over £/ e 's ( a positive v e v e asymmetry), which could lead to a reduced neutron-to-proton 
ratio and, hence, a lower Y. 

Previous bounds on these mixing channels have been calculated assuming a small r\ 
(« 3 x 10~ 10 ) because an accurate measurement of the primordial D/H was lacking. It is 
now, however, established that the primordial D/H lies at the low end of the range previously 
thought: D/H^ 3.4 ± 0.3 x KT 5 pf. This low value infers a higher 77 « 5.2 ± 0.3 x 10~ 10 , 



which predicts a slightly higher primordial He abundance Y for the standard picture N u = 3. 



The upper limit on N v based on the measured Y should then become tighter [|T(J . 

However, complicating this issue is the debate over the systematic uncertainties in the 
observationally-inferred Y . Based on 62 low metallicity HII regions, Olive et al. |17j esti- 



mated Y = 0.234 ± 0.002 (stat.) ± 0.005 (sys.); or Y = 0.237 ± 0.003 (stat.) ± 0.005 (sys.), 
when the controversial object I Zw 18 (also the most metal- deficient object of the sample) is 
ignored. Based on another sample which partially overlaps with that of Olive et al, Izotov 
and Thuan |T8| have claimed Y = 0.244 ± 0.002, partly because they deduced a higher Y 



for I Zw 18. Disagreement aside, the two values do not seriously contradict each other given 
the systematic uncertainty ~ 0.005 of the estimates. With the systematic uncertainty in 
mind, we can still set a conservative upper limit Y < 0.25 for the purpose of constraining 
the neutrino oscillation physics at the BBN epoch. A higher Y would not only be inconsis- 
tent with both estimates, it could also possibly be inconsistent with the morphology of the 
Horizontal Branch stars in globular clusters . 



Thus, given D/H^ 3.4 ± 0.4 x 10 5 and Y < 0.25, we can set a rather conservative 
bound on N v at the 95% C.L., 

N v < 3.3. (8) 

This bound can be used to constrain the active-sterile neutrino transformations in the 
BBN epoch, but it is more precise and convenient to calculate directly the helium yield of 



neutrino oscillations, because it is ultimately the helium mass fraction that is compared to 
observations. We took rj = 4.6 x 10~ 10 and recalculated the 4 He yield in the presence of 
active-sterile neutrino mixing. (Our calculations are similar to those in Shi, Schramm and 
Fields ]7|.) Figure 1 shows the result. A parametrization of the bound on the <-> v s 
mixing can be obtained: 

\5m 2 \ sin 4 29 < 1(T 5 (1(T 7 ) eV 2 for 5m 2 > (<)0. (9) 

The <-> u a mixing solution to the SuperK atmospheric neutrino data is clearly ruled out 
by BBN if lepton asymmetry is neglig ible (L < KT 7 in this case). To put this result in 
perspective, our results suggest that a v s mixing solution to the SuperK data in this 

scenario would imply N v = 3.9. 

IV. BBN BOUND IN THE PRESENCE OF AN AMPLIFED LEPTON 

ASYMMETRY 

A resonant v a <-> v s transformation can amplify the lepton asymmetry in the v a sector 
to a level 

T __ 3(rc„ a -n Da ) |£m 2 /eV| 

K a = ± ^w?~' ( 10 ) 

where Tq = T/MeV @. Generation of this lepton number does not necessarily force v s into 

chemical equilibrium if the mixing is small enough to satisfy the constraints in Fig. 1 ||. 

To calculate the amplification of the lepton asymmetry by v a <-> v s mixing, we take the 

initial lepton number to be negligibly small. Note that in this limit L Va = 3(P Z — P Z )/1Q. 

From Eq. and its anti- neutrino counterpart, L Va satisfies the approximate equation 

dL Ua DV[ f 3V f3(2L Ua + L )P Z 

dt ~ V 2 + [V - (5{2L Va + L )Y 1 4{V 2 + [V + (3(2L Ua + L )] 2 }L 

which is valid in the adiabatic limit, |V| > |dV/dt|/|V|. In Eq. (01]) we use V = (V z + V z )/2, 
and (3 = (V z — V z )/2(2L Ua + L ) w 0.35GfT 3 . We have also assumed P z ~ 1, because the 
v a *^ v s mixing would violate the BBN bound in Fig. 1 if P z <^ 1. 
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Eq. (|Tl|) is a damping equation for 2L Ua +L , unless Vq > 0. Here Vq > if 5m 2 < (i.e., 
a heavier z/ Q ) and if the temperature T drops below the resonant temperature for v a <-> za, 
(see the Appendix), 

<5m 2 



22 



1/6 

MeV for a = /i, r. (12) 



leV 2 , 

For T > T res , 2L Ua + L will be driven toward zero. In this case approaches zero, 
unless the initial asymmetry Lq (whose definition may be expanded here to also include the 
neutrino asymmetry generated by other active-sterile neutrino mixings) is too large for the 
damping to be efficient. Efficient damping occurs when the generated lepton asymmetry in 
other neutrino species satisfies 

\L\ < 10- 4 |5m 2 sin2^| 1 / 2 , (13) 

so that (dL Va /dt)/ L Va < H, where H is the Hubble expansion rate. 

Eq. (Q) and its anti-neutrino counterpart are non-linear equations. It has two attractors 
which give the asymptotic values of Eq. (|T0|). The asymmetry L Ua calculated through nu- 
merical integrations of these equations exhibits chaotic behavior for T just below T rcs 0. In 
this phase L Va oscillates chaotically between the two attractors. For sufficiently low T, L Va 
approaches asymptotically one of the two attractors. It is in this asymptotic phase when 
the adiabatic condition is satisfied and Eq. (|TI]) applies. 

It is conceivable that a v T <-» v s i mixing satisfying 

r ,2 1/6 

sin 2 2^>10- n . (14) 



eV 2 



and the the BBN bound Eq. ([J), can generate \L Ut \ ^> 10 7 P,|20|, which in turn suppresses 



the bad BBN effects of *-> v s transformation with 5m 2 ~ 10 -3 eV 2 and sin 2 29 ~ 1, and 
so relaxes the BBN bound on the <-> v s mixing. (Here u s > and v s could be in principle 
the same species but are distinguished here for clarity.) But as we will demonstrate below, 
while this is not impossible, there are severe constraints on the required v T «-> v s > mixing 
parameters for this loophole to be realized. 
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To distinguish between the v T <-> v s i transformation and the <-> z/ s transformation, all 
quantities applying to the former process are marked with an prime. To suppress the trans- 
formation due to a near maximal <-> v s mixing, the lepton asymmetry has to be amplified 
early by the v T <-> u s > transformation channel before any significant <-> z/ s transformation 
can occur. This requires the v s production rate at the v T «-> zv resonance temperature T r ' es 
to satisfy 

„ (5m 2 sin 20 \ 2 rr . . 

D HH <ff ' (15) 

with H = 5.5T 2 /m p \. Here m p \ ~ 1.22 x 10 28 eV is the Planck mass. Employing T r ' es and 



Vj we obtain from Eq. fll5|) : 



|<W 2 | > 10 1 5m 2 sin 20 1 4/3 ~ lO^eV 2 . (16) 
Once the lepton asymmetry is amplified, the <-> z/ s mixing will have 

() 771 ^ $771 ^ 

^ = -W Sin2 ^' ^ = -^ cos20 + T f ± ^~ Vji^L^ + L.J, (17) 

where "+" ("— ") applies to the z/^ <-» za, (z/^ <-» p a ) transformation channel. If L Vt is positive 
(negative), the <-» v s (z? M <->• za,) channel will experience a resonance, regardless of the sign 
of 5m 2 , at a neutrino energy 

R _ P\L Ut \ _ [Luj 
res ~ V£/E ~ 63G F T 1 j 

(assuming initially L„ = 0). The resonance moves from lower energies to higher energies 
in the mu neutrino spectrum as \L Vt \ grows. At the same time, the u s [y^ <-> v s ) 

transformation channel is suppressed by L Vt as a consequence of having a even larger \V Z \. 
The net result is a newly generated L v with a sign opposite to that of L Vt . If this L u 
becomes as large as ~ —L Ut /2, the «-> v s mixing will be unhindered and will eventually 



bring v s into equilibrium. Based on this consideration, Foot and Volkas pO| , pT| argued for a 
requirement on the v T <-> v s mixing 



\5m' 2 \ > 20 eV 2 . (19) 



Here we show that this requirement is too weak. We find that a more stringent requirement 
can be obtained based on the same consideration at the beginning stage of the L Ut growth, 
when both L Vt and the mu neutrino resonance energy E ms are very small. For definitive- 
ness, we assume that L Ut > so that both the v T <-> v s i transformation and the <-> v s 
transformation encounter resonances. 

Each of these two resonances occurs in a finite neutrino energy bin (resonance width). 
Outside the resonance energy bins the effective mixing between active neutrinos and sterile 
neutrinos quickly diminishes and becomes negligible. This feature enables a simple semi- 
analytical approach to accurately calculate the growth of lepton number asymmetries by 
tracking only neutrinos in resonance regions. The details of this semi-analytical approach 
are presented in the Appendix. 

The constraint on 5m' 2 results from the requirement that \L„ T \ be much larger than \L V \ 
at any temperature as the <-> v s resonance sweeps through the entire neutrino energy 
spectrum. In other words, we must have 

^rcs 



/(e res )5e res R 



< \ (V + 2L„„) (20) 



de rcs /dt 

for any e rcs - In the equation, /(e) = [2/3£(3)] [e 2 /(e e + 1)] is the neutrino distribution 
function. <5e res is the energy width of the resonance. f(e Tes )Se ies is therefore the fraction of 
mu neutrinos undergoing resonance. R is the resonant transition rate. <5e r cs/|de rC s/dt| is the 
duration of the resonance at e res . The energy width of the resonance depends on whether 
the resonant transition is collision-dominated (D > V x ), or oscillation-dominated (D < V x ): 

2\Dde xcs /dV z \ if D > V x 



5e rns ~ < 



(21) 

2\V x de Tes /dV z \ ifD<V x 



Likewise, the resonant transition rate depends on these parameters: 

V x 2 /D if D > V x 



R 



(22) 

V r ifD< V T 



We only consider the case where e res 1, because this is when L Vt is in its initial stage of 
growth [L Vt <^ 10 -7 ) and is most easily matched by a competing L u ^. In this case, T T r ' es 
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(the temperature at which L Vt growth starts), and /(e res ) ~ e 2 cs /1.8. We can further rewrite 
|e res | = ife res |dlne r es/dlnT|. Then with Eq. (|I2"D and the expression for H, Eq. (|20| ) becomes 

5m' 2 /leV 2 \ U/6 > 2 x 10 4 e r ~i (dlne^/dlnTf 1 |<5m 2 /10- 3 eV 2 | 2 if D > V x ; (23a) 
5m ,2 /leV 2 | 17/6 > 10 3 e r -J (dlne^/dlnTf 1 jcW/lO- 3 eV 2 | 3 if D < V x . (23b) 

In the initial stage of rapid L Vt growth, |dlne res /dlnT| can be related to the growth rate 
of L Vt by |dlne res /dlnT| ps |dlnLj, T /dlnT — 2| |dlnLj, T /dlnT| (with L„ safely ignored). 
Our semi-analytical calculations show that in the initial exponential stage of L Vt growth 
when L Vt is <C 10 -7 , |dlnLj, T /dlnT| is well approximated by 

d In L„ 



dlnT 



6 x 10 6 sin2fl', 



(24) 



indendent of 5m' 2 . Of course, sin 26' here must satisfy the BBN bound Eq. (|9|). 

Eq. (23) and ( 24]) show that the most stringent requirement on m 2 r — m 2 ; does indeed 
come not from <-» z/ s resonances at e res ~ 3, but from resonances centered at the smallest 
possible e res as long as the or transition probability in that resonance energy bin is 
<C 1. This condition, expressed as 

5e rc , 



R 



de res /dt 



< 0.1, 



(25) 



can be rewritten as 



c > 

-res 



8m' 2 


-1/2 


d In e res 


-1/3 


5m 2 


2/3 


leV 2 




dlnT 




10" 3 eV 2 





(26) 



regardless of the value of D/V x . It can be shown that e res is in the oscillation-dominated 
regime if 



Hres £ 0.25 



5m 2 


1/2 


dm 


-1/2 


10" 3 eV 2 




leV 2 





(27) 



Therefore, the most stringent requirement on 5m' 2 comes from the oscillation-dominated 
regime for sin 2 26' > 10~ 8 (while |dlne res /dlnT| > 10 3 ), and from collision-dominated 
regime for sin 2 26' < 10~ 8 (while |dlne res /dlnT| < 10 3 ). 
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Combining Eq. (23), Eq. fl24|) and Eq. ( p6|) yields a requirement on the mass-squared- 
difference necessary to effect suppression of *-> v s transformation at the Super Kamiokande 
level: 

200 ( \ml - ml |/10- 3 eV 2N ) 3/4 eV 2 for sin 2 26' > lO" 8 



<-ml f >\ v -< / ~ (2g) 

(sin 20') /2 (|^ M - m 2 J/10- 3 eV 2 ) eV 2 for sin 2 29' < 10~ 8 

Since |<5m' 2 | = m 2 T — m 2 Vs , Eq. ( |28| ) implies that we must have 

m, T >15eV (29) 

to successfully suppress the <-» v s transformation in BBN. 

Tau neutrinos this massive would contribute a fraction of the critical density Q u ~ 0.5h^Q 
(where h^o is the Hubble constant in the units of 50 km/sec/Mpc) in the form of Hot Dark 
Matter (HDM) ||22|| . When normalized to yield the observed structure today, a matter- 
dominated flat (Q m = 1) model universe with Q u > 0.3 (the remainder being mostly Cold 
Dark Matter) yields too few Damped Lyman-a systems (protogalaxies) at a redshift z > 3 
to accommodate observations |2^J2^] . Due to the free-streaming of neutrinos, the Hot Dark 
Matter component reduces the density fluctuation amplitudes at the galaxy mass scale and 
causes galaxy-sized structures to form too late [z < 1). Even an Q u = 0.2 flat, matter- 



dominated model may still be in disagreement with observations |[25|| . Therefore, any light 



stable and weakly-interacting neutrino more massive than m v > 5/i| eV will have trouble 
with structure formation considerations. Models with these neutrinos cannot be rescued by 
pushing hso much higher, because both the age of the universe and the observed structure 
today require h^o pa 1 [26|. If as implied by the cosmic deceleration measurements [^7|, the 



total matter density in our universe is less than the critical density, the incompatibility of 
structure formation at high z and Q u > 0.2 will become even worse - the HDM would then 
comprise a larger fraction of the matter density and so reduce the density fluctuations at 
small scales even more effectively. 

To simultaneously have m„ T > 15 eV and successful structure formation at high redshifts, 
the massive tau neutrinos would have to decay before making an imprint on the cosmic den- 
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sity fluctuation spectrum. Since sub- horizon-sized density fluctuations (i.e., those affected 
by free streaming of neutrinos) were frozen before the epoch of matter-radiation equality 
and grew only after that, the density fluctuations residing in other Dark Matter components 
would not be affected by the density fluctuations in the tau neutrino spatial distribution so 
long as these z/ r 's decayed before the matter-radiation equality epoch. The lifetime of u T 
must then be 

r„ T < 10 3 year. (30) 
The decay cannot be radiative, as it would violate the bound on neutrino electromagnetic 



dipole moments |2 



For those u T «-> u s mixings with 10 2 eV 2 < |<5m' 2 | < 10 1 eV 2 , the subsequent 



v., <-> 



u s resonance occurs during the chaotic phase of the u T <-> u s > evolution. Our numerical 
calculations show that in this case the two mixing systems are coupled and continue to 
be chaotic. The u^ <-> u s transformation will in fact soon dominate the chaotic oscillation 
process because this transformation channel has a larger \V X \ than the BBN-constrained 
u T <-> u s i transformation channel. Therefore, for 10~ 2 eV 2 < |<5m /2 | < 10 _1 eV 2 , the u^ <-> 
u s transformation cannot be suppressed by the presence of the u T «-> u s > transformation 
channel. 

The u T <-> u s i system is also bounded by BBN at large 5m' 2 . If a large (L Ur > 0.1) 
asymmetry is generated by the u T <-> u s i transformation before T ~ 5 MeV, the produced v s i 
or u s i and the re-thermalized u T u T will increase N v by 0.5, violating the BBN bound. This 
imposes an upper bound on \5m' 2 \ |T0| . p9| : 



\Sm' 2 \ < 10 4 eV 2 , (31) 

implying m Ur < 100 eV if u s > is much lighter. 

In figure 2 we plot the conditions that u T and the u T ^ u s > mixing must satisfy in order to 
alleviate the BBN constraint on the maximal vacuum u^ «-> u s mixing that may be implied 
by the SuperK data. 
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We finally very briefly comment on another scenario in which a positive lepton asym- 
metry generated by a v T <-> v s i mixing may be partially converted into the v e v e sector by 
a matter-enhanced v T <-» u e transformation. This positive u e chemical potential can reduce 
the effective N v and may make some room for extra neutrino energy density A previous 
paper has argued that the reduction in N v can be as large as AN U ~ 0.5 [ID]. This will not 
accommodate for an extra AN V = 0.9 produced by the <-> v s mixing in question, given 
the 4 He-derived bound N„ < 3.3. We also note that even the 0.5 reduction in N v may be 
overestimated [p9|. 



V. SUMMARY 

We have shown that based on the observationally-inferred primordial abundance of 4 He 
(Y < 0.25) and deuterium (D/H« 3.4 ± 0.3 x 10 -5 ), Big Bang Nucleosynthesis yields a 
stringent bound on the effective number of neutrino species (energy densities) during the 
BBN epoch, N v < 3.3. This bound can be employed to constrain active-sterile neutrino mix- 
ings, as plotted in Fig. 1. With the preassumption that all three active neutrinos are light 
compared to 1 MeV, it rules out the ^ v s mixing explanation of the Super Kamiokande 
atmospheric neutrino data. The only way to circumvent this bound is if there is a simulta- 
neous mixing between v T and a lighter sterile neutrino v s > with the following properties: 

1. A mixing mass-squared difference 200 eV 2 < m 2 ^ — m 2 / < 10 4 eV 2 and a non-radiative 
v T decay lifetime t Vt < 10 3 year. The mass-squared difference implies a u T mass 
between ~ 15 eV and ~ 100 eV if u s i is much lighter than v T . 

2. [(ml T -m 2 s ,)/eV 2 ] 1 / 6 sin 2 2^ > lO" 11 and [(m 2 T - ml si )/eV 2 ] sin 4 26 < 10" 7 , where 9 
is the vacuum mixing angle between v T and v s <. 

We conclude that the <-> v T vacuum oscillation channel with maximal or near maximal 
mixing angle provides a more natural explanation for the Super Kamiokande data. 
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VI. APPENDIX 



Our semi-analytical calculations track the movement (in energy space) of resonances in 
the neutrino energy spectrum and the resonant conversion rate of neutrinos in the resonance 
regions. For v a «-> u s mixing with vacuum mixing parameters 5m? = m? Us — m? Ua < and 
sin 2 29, both the conversion of v a and v a are calculated. If we denote as the resonant 
neutrino energy (divided by the temperature) for v a ( "+" ) and v a ( "— " ) , The rate of resonant 
conversion, when normalized by the total u a or u a number density, is 



de (±) 



(±) 



V.' = rf[^\^^«. (32) 



This equation is essentially the l.h.s. of Eq. ( P0f) except for a factor tt which comes from a 
detailed integration of transition probability over the entire resonance region. The rate of 
change of L Va is therefore 

^ = 2 fo<+) _ K^) . (33) 
dt 8 V > K ' 

When L Va is small (i.e., \V£\ *C |V^|), the effect of matter-antimatter asymmetry on the 
v a *^ v s and v a <-> v s oscillation is small. Both oscillation systems will encounter resonances, 
at an energy 



e 



x ± 1 35G F T 3 (2L 1 , Q + L ) 



\5m 2 \/2T 



(34) 



(see Eq. |§ to [|]). Eq. ([33|) is approximately 

* i^t GfT3 sin2 29{L ^ + Lo/2) ' (35) 

with d/(e)/de evaluated at e = |5m 2 /44G 2 7 T 6 | 1 / 2 . This is apparently a damping equation 
for (L Ua + L /2) if d/(e)/de < (when e < 2.217), and a growth equation for (L Va + L /2) 
if d/(e)/de > (when e > 2.217). Therefore, (L Ua + Lq/2) is always damped to at 
temperatures T > T r ' es , when is small; and (L Va + Lq/2) grows for T < T r ' es when is 



large enough. 
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In Figure 3 we plot the result of our numerical calculation based on Eq. (|32| ) and (|3^) for 
v T *-> v s i mixing, assuming m 2 T — m 2 f = 50 eV 2 and sin 2 29' = 10~ 4 . It is for the most part 



similar to the thick solid line in Figure 1 of Foot |2l], which employs the same parameters. 
There are minor differences that are readily identifiable: (1) our result tracks T -4 more 
closely in the "power-law growth" epoch; (2) L Vr in our results does not switch sign at the 
initial point of growth. The sign difference is not surprising because of the chaotic character 
of the growth, which introduces a sign ambiguity to the problem |§. Similar calculations 
for other choices of m 2 T — m 2 / and sin 2 29' = 10~ 4 show that the temperature at which the 
L Vt growth starts is approximately given by Eq. (|12|). 

As an illustration, also plotted in Figure 3, are the \L Vr \ required for the «-> v s 
or <-> v s resonance to occur at or energies e res = p^l/T = 0.01, 1, 10. (The 
parameters for the <-> v s mixing are 8m? = 10 -3 eV 2 and sin 29 = 1.) It can be seen 
from the figure that the lower energy component of the or neutrinos encounters the 
resonance first when \L Vr \ is very small, and the resonance region moves through the or 
Up spectrum to higher neutrino energies as \L Vr \ becomes larger. Essentially all z/ M or z/ M 
encounter resonances at T T r ' es . 

In Figure 4 we show |dlnLj, T /dlnT| as a function of the v T <-> v s > vacuum mixing 
parameters, m 2 — m 2 s and sin 2 29', in the initial exponential stage of L Ut growth when 
L Ut is 10 -7 . This figure shows that |dlnLj, T /dhiT| is approximately a linear function 
of sin 29', and is insensitive to m 2 T — m 2 s . An analytical fit to this numerical result yields 
Eq. (m. 
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Figure Captions: 

Figure 1. 4 He yields of the (orz/ T ) <-> v s mixing (for r] = 4.7 x 10~ 10 ). Mixings that yield 
Y > 0.25, including the region inside the dashed line that is compatible with the Super 
Kamiokande data, are ruled out by BBN. 

Figure 2. The parameter space of v T <-> v s i mixing and the additional requirements on 
v T that may alleviate the BBN constraint on the <-> v s mixing parameters which can 
accomodate the Super Kamiokande data. 

Figure 3. The growth of the tau neutrino asymmetry as a result of the tau neutrino-sterile 
neutrino mixing, assuming m 2 Vr — m 2 i =50 eV 2 and sin 2 29' = 10~ 4 . The intersections 
between the growth curve for L Ut and the dashed lines indicate when resonances occur for 
(if L Vt > 0) or Dp (if L Ut < 0) neutrinos with momentum p. 

Figure 4. The initial rate for L Ut growth, dlnL/dlnT, as a function of the vacuum tau 
neutrino-sterile neutrino mixing parameters. 
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